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ABSTRACT 

The unveiled main-sequence splitting in uj Centauri as well as NGC 2808 suggests that matter 
highly-enriched in He (in terms of its mass fraction Y 0.4) was produced and made the color of 
some main-sequence stars bluer in these globular clusters (GCs). The potential production site for the 
He-rich matter is generally considered to be massive AGB stars that experience the second dredge-up. 
However, it is found that massive AGB stars provide the matter with Y ^ 0.35 at most, while the 
observed blue-shift requires the presence of y ~ 0.4 matter. Here, we show that extra mixing, which 
operates in the red giant phase of stars less massive than ~ 2 Mq, could be a mechanism that enhances 
He content in their envelopes up to y ~ 0.4. The extra mixing is supposed to be induced by red giant 
encounters with other stars in a coUisional system like GCs. The Y ~ 0.4 matter released in the AGB 
phase has alternative fates to (i) escape from a GC or (ii) be captured by kinematically cool stars 
through encounters. The AGB ejecta in uj Cen, which follows the latter case, can supply sufficient 
He to cause the observed blue-shift. Simultaneously, this scheme generates the extreme horizontal 
branch, as observed in lu Cen in response to the higher mass loss rates, which is also caused by stellar 
encounters. 

Subject headings: globular clusters: general — globular clusters: individual {uj Centauri) — stars: 
AGB and post-AGB — stars: evolution — stars: horizontal-branch 



1. INTRODUCTION 

The discovery of a split in the main sequence (MS) of 
UJ Cen (jBedin et al.ll200l ) has opened a new window in 
the study of Galactic globular clusters (GCs). The ob- 
served fact that stars of the blue MS (bMS) in uj Cen 
exhibit slightly stronger absorption of iron l ines on aver- 
age , in comparison with the red MS (rMS) (jPiotto et al.l 
|200E| ) which strongly implies that the origin of bMS is 
attributable to the enhancement of He inside bMS stars. 
Subsequently, the recent HST observation has provided 
us with the second sample of MS splitting, NGC 2808 
(jPiotto et all 120071 ). This GC essentially exhibits no 
dispersion of [Fe/H] and thus the argument of He en- 
hancement in bMS stars has been reinforced. A com- 
parison between the observed color-magnitude diagrams 
(CMDs) and synthetic population models suggests that 
bMS He abundance is enhanced t o Y ~ 0.4 for both of 
two GCs (iNorris 2 004: Lcc ct all 120051: iD'Antona et all 
120051 iPiotto et al.l[2007t ). iTsujimoto et all (|2007f ) have 
argued that if the surface of a star on the MS is polluted 
with the Y = 0.4 matter with the mass of ~ O.IA/q, the 
star can move to a position on the observed bMS in the 
CMD. This relaxes a severe demand on the amount of He 
possibly supplied from asymptotic giant branch (AGB) 
stars, in contrast to the idea that the bMS stars are born 
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from pure AGB ejecta c onsist ing of y 0.4. Similarly, 
iNewsham fc Terndrup I ()2007[ ) have insisted that surface 
pollution may explain the high He content of bMS stars. 

It should be stressed that Y ~ 0.4 matter is necessary 
to realize the bMS stars in these GCs. Since massive 
AGB stars can enhance the He abundance in their en- 
velopes through the second dredge-up in the early AGB 
phase, they have been con sidered a possible pro duction 
site for the Y ~ 0.4 matter (jD'Antona et al.ll2005[ ). How- 
ever, previous studies indicated that the resultant abun- 
dance of He in the envel ope of any AGB star does not 
excee d Y ~ 0.35 (see e.g.. lvan den Hoek &: GroenewegenI 
Il997t l. Therefore, such a He yield from massive AGB 
stars may not be sufficient to split the MS as observed. 
In the subsequent section, we will discuss this issue and 
conclude that it is highly implausible for a massive AGB 
star to enhance He abundance up to y ^ 0.4 in its enve- 
lope. 

Then, we should search for an alternative mechanism 
that enhances the He abundance more efficiently than 
the second dredgc-up. Wc believe that this process, 
which enhances He abundance, should be accompanied 
by some other elemental signatures. This is reminiscent 
of the abundance anomaly of red giants in GCs. The 
abundance anomalies observed for red giants, such as 
the abundance variations of CNO elements and the 
0-Na anticorrclation, arc common att ributes associ- 
ated with GCs including uj Cen (see iGratton et al.l 
|2004[ ). One of the proposed mechanisms to explain 
these anomalies is a deep mixi ng, so-called extra 
mixing, on the red giant s tage (ISweigart fc Mengell 
19791: iLanger et all Il993l: ICharbonnel at al.l Il998l: 



Fuiimoto et al.lll9'9^lAikawa et al.ll2001l : IChaname at al.l 



20051 ISuda fc Fuiimotoll2006[ ). The extra mixing is as- 



sumed to be caused by some kind of rotation-induced 
mixing, driven by shear instability around the border 
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of the He core (see iFuiimotol [l988l: IZahnl [l99l) . Since 
such a deep mixing naturally draws He from the core, 
the extra mixing is a promising mechanism of producing 
He-rich matter. As a mechanism of mixing, most of the 
previous works assumed a continuous mixing caused 
by the meridional circulation, while Fujimoto and his 
coworkers insisted that a flash of H brought into a 
degenerate core by the rotation-induced mixing drives 
the intermittent mixing. 

In this study, we propose that extra mixing, which oc- 
curs in the red giant phase of a star, enhances He abun- 
dance in the envelope and ejects He-rich matter when the 
star evolves to a white dwarf through mass loss. In this 
scenario, the driving force of the extra mixing is gen- 
erated through interactions between the red giant and 
dwarfs in the central region of to Cen. Our detailed cal- 
culations reveal that stars with the masses less than ~2 
Mq can enhance He abundance up to F ~ 0.4, provided 
that the extra mixing operates in the red giant phase. 

Based on our scenario, we discuss other conspicuous 
characteristics of to Cen, the extremely blue horizontal 
branch (HB), which is observed in some other Galactic 
GCs. Although factors that determine HB morphology 
is still an open question and is known as the second pa- 
rameter problem, the most effective parameter that af- 
fects the location of a HB star in the CMD is the mass 
Menv of the envelope. A major factor to determine M^nv 
is the mass loss during the red giant phase. We predict 
that encounters induce a high mass loss rate through a 
gain of angular momentum and promote the presence of 
an extremely blue HB population. In the end, two mys- 
teries in uj Cen of MS splitting and HB morphology arc 
discussed in a unified framework of encounters of dwarfs 
with AGB ejecta and those of red giants with dwarfs. 

2. MASSIVE AGB STARS AS A SOURCE OF HE-RICH 
MATTER 

In this section, we estimate He abundance in the enve- 
lope of a massive AGB star after the second dredge-up 
and show that it is enhanced as high asY ~ 0.35 but can- 
not significantly exceed this value. The He enhancement 
Ayonv in the envelope due to the second dredge-up is de- 
termined by (1) the envelope mass, (2) the He core mass 
Mcore just before the second dredge-up, and (3) the max- 
imum depth of the convective envelope during the second 
dredge-up whose mass coordinate is given by Mmix.basc- 
For a massive AGB star, the envelope mass is not reduced 
significantly by the mass loss until the beginning of ther- 
mal pulses. The surface He abundance becomes smaller 
for a more massive envelope since the matter dredged-up 
from the He core is diluted in the whole convective en- 
velope. Assuming a simple stellar structure consisting of 
the He core and the H-rich envelope, the He enhancement 
Ayonv in the envelope is approximated by the following 
formula; 

A (-^-^coic -^^mix, base) ( 1 ^nv) 

-'^-'init -^■'-'mix, base 

where Mian denotes the initial mass of a star. The core 
mass Mcore is calculated with the stellar evolution pro- 
gram for a set of given initial values of mass. The initial 
value of Ycnv is assumed to be 0.24 in the rest of this sec- 
tion. Using this formula, we can estimate the amount 
of the dredged-up matter during the AGB mass loss 



phase to attain a given enhancement of He. For example, 
Mmix,base ~ 0.32 Mq is required to attain Al^nv — 0.15 
by the second dredge-up for a star with Minit = 7 Mq 
(Mcore — 1-68 Mq) if thc envelope does not significantly 
lose its mass. Of course, it is impossible to dredge-up He 
from such a deep interior since thc C-0 core of this par- 
ticular star becomes as massive as 0.9 Mq before the 
second dredge-up. Therefore, this star cannot enhance 
the He content in the envelope to Yenv = 0.4. 

After thc second dredge-up, thc mass of He in the He- 
burning shell is reduced from inside due to the growth of 
the C-0 core. Thinning of the He-burning shell inhibits 
further enhancement of the He content in the envelope by 
the subsequent third dredge-up. For example, in a 7 Mq 
star. He shell flashes convert some He into C and dredge 
up the He layer with thc mass of 0.1 Mq (Y ~ 0.8) to the 
envelope. The He enhancement for a 7 Mq star (6 Mq 
envelope) is estimated to be only AY^nv ^ 0.013. As a 
result of the above formula as well as previous studies on 
stellar evolution ;Hc content in the envelope of a massive 
AGB star with Mi^i > 3 Mq can attain Y^nv ^ 0.36 
at most and never exceeds 0.4. Indeed, no established 
models have been proposed to attain Y^^v ^ 0.4 in the 
envelope of a massi ve AGB star. As anot her problem, it 
has been discussed ([Lattanzio et al.ll2004l ) that the yield 
of CNO and light p-capture elements in massive AGB 
stars cannot explain the abundance anomalies in GCs. 

3. STELLAR MODELS WITH EXTRA MIXING 

Models in this work are com puted with a stellar evolu- 
tion program based on llben e t al. (1993), wh ich includes 
updated ph ysics and improved opacity tables (jSuda et al.l 
12004 120071 . Suda & Fujimoto 2007, in preparation). To 
construct isochrones for MS, we evolve stars through 13 
Gyr starting from the zero-age MS for various initial 
masses at an interval of 0.01 Mq. Each evolutionary 
point is interpolated from the lower MS to subgiant. On 
thc other hand, thc isochrones for the red giant branch 
(RGB) are substituted by the evolutionary tracks of a 
star with the appropriate mass, since the isochrones on 
the RGB coincide with thc evolutionary tracks for stars 
with thc same initial masses and metallicities. In this 
RGB phase, we assume that the mixing event occurs 
through encounters in a GC and the matter in the core 
is intermittently dredged up to the envelope and mixed. 

The total amount of He in the envelope is determined 
by the mixing rate, which is the amount SM^i^ of mixed 
He per event divided by thc interval Tmix between the 
mixing events. The duration of the RGB phase should 
be much longer than the time for a star to cross the 
cluster core D/v ~ 1 pc/lOkm s""'^ 10^ yr, where D 
denotes the size of the cluster core and v the velocity 
of a star. While a red giant passes the cluster core, the 
star often encounters other stars and its envelope ac- 
quires some amount of the orbital angular momentum. 
Thc resultant excess rotation in thc envelope causes mix- 
ing through shear instabilities at the boundary between 
thc core and envelope. There is a critical value for the 
mixing rate (JMmix defined by the dredged-up mass per 
mixing event; (SM„iix/Tmix = (LH)/(£^H^env), where 
denotes the energy per unit mass released by H-burning, 
(in) the average luminosity of the star, and X^nv is the 
H abundance in the envelope. The righthand member 
increases as helium mixing proceeds. In our models, the 
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values of JMmix arc chosen to be near the obtained crit- 
ical mixing rate. If (S A/mix is too large relative to the 
critical rate, the He core cannot grow. As He mixing 
proceeds, however, it grows since the critical rate is a de- 
creasing function of Xcnv The border of the occurrence 
of He mixing is assumed to be Mini = 2.1Mq. A star 
less massive than this has a degenerate He core. Once 
such a star on the RGB encounters another star, the 
angular momentum in the envelope increases to induce 
shear between the core and envelope. Then a fraction of 
H in the envelope is repeatedly mixed into the core and 
undergoes the off-center flash. T his off-center H flash i s 
believed to induce extra mixing (iFujimoto et al.lll999f ). 
In more massive stars, the mixed H would not undergo 
a flash because electrons in the core are not degenerate. 



4. RESULTS 



4.1. 



Y as a function of stellar mass 

The evolutions of He core mass and content yonv in 
the envelope during the RGB for a star with Afini = 
1.5Mq are shown as an example in Figure 1. Basically, a 
larger JAfmix leads to greater He enhancement. For the 
case with 5Mmix= 0.0015 Mq, the core that was once 
reducing its mass begins to grow after the flrst several 
mixing events due to the increasing l^nv Then, it grows 
rapidly as the H-burning rate increases with the core 
mass. The resultant Y^nv at the tip of RGB becomes 
greater than 0.4. This, together with other low-mass 
stellar models, reveal that stars with A^ini < SM© can 
achieve Y^nv > 0.4 at the tip of the RGB as a result of 
the extra mixing. These He-rich matters are ejected at 
the final stage of AGB. Figure 2 shows the He abundance 
l^nv finally attained in the envelope as a function of the 
initial mass Afini of stars in the range of O.SM© < Afini < 
8Mq , together with the average yield from massive AGB 
stars with Afini > 3Mq by other studies. Note that Y ~ 
0.4 matters are produced from low-mass stars only in 
GCs where red giants encounter dwarfs. Additionally, it 
should be stressed that the He-rich matter is expected 
to show C depletion and N enhancement as a result of 
partial CN cycles. For the parameter range surveyed by 

this work, we obtain [C/Fe] = -0.2 0.5 and [N/Fe] 

= 0.8 — 0.9. These results are marginally compatible with 

the observed values [C/Fe] 0.1-0 and [N/Fe] > 1 

for bMS stars (Piotto et al. 2005) and subgiants with 
<[Fe/H]>=-1.37 (jVillanova et al. 2007). 

On the other hand, stars with A/ini > 2.1Afo enhance 
the He abundance in their envelopes through the process 
to deepen a convective envelope by the second dredge- 
up. Furthermore, stars with a mass of ~3.0 Mq stay on 
the RGB only for ~ 10^ yr. Therefore, the chance of 
an encounter is so limited that such stars are unlikely to 
undergo extra mixing, even if the extra mixing did not 
need a degenerate He core. The Vonv of these massive 
AGB stars increases with increasing stellar masses and 
a A/ini = SMq star giv es yenv 0.35 as the maximum 
yield. As presented bv iHerwigl (|2004[ ). a model adopt- 
ing overshooting leads to efficient enhancement of l^nv 
after the second dredge-up. At the same time, however, 
it accelerates the growth of the AGB core. As a result, 
the lower mass limit of a star that undergoes central C- 
burning decreases to A/ini ^ TMq from ~ 8 A/©. There- 
fore, even with overshooting, it is unlikely to realize l^nv 
~ 0.4 in massive AGB stars. 
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Fig. 1. — Time evolution of the He core mass A/corc and the He 
abundance Vonv in the envelope during the RGB phase for a 1.5 
Mq star without (dotted line) and with extra mixing. In the extra 
mixing models, iSA/mix = 0.0014 Mq (dashed line) and 0.0015 Mq 
(solid line) are assigned. The value of r^i-^ is set to be 3 X 10^ yr. 
The calculations are performed until the tip of the RGB. 
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Fig. 2. — He abundance ycnv finally attained in the envelope 
as a function of the initial stellar mass, Afinj. For Afini < 2.1AfQ, 
extra mixing is assumed to occur during the RGB. The interval of 
mixing is assumed to be T^ix = 3 X 10^ yr. The values of the mixed 
mass per event are set to be &M^^^= 0.0015, 0.0025, 0.0035 Mq , 
for models with Afinit = 0.8, 1.5, 2.0 Mq, respectiv ely. Our results 
for Afjni > 'AMq are compa r ed wi th other studies bv lHerwig I II2004I) 
and I Ventura fc D'Antona I II2005I) . 



4.2. Surface Pollution and Main-Sequence Splitting 

As discussed in iTsuiimoto et all l)2007f ). a 0.1 A/© ac- 
cretion of y = 0.4 matter by each star is sufficient to 
reproduce a bifurcation of the MS as observed in uj Gen. 
Thus, the total He mass supposed to be accreted by the 
bMS stars is estimated to be ~ 0.012A/(^, where A/^^ de- 
notes the stellar mass of uj Gen, if the mass fraction of 
the bMS stars is assumed to be 0.25. On the other hand, 
the amount of He ejected from low mass AGB stars of 
all generations of stars is calculated to be 0.02A/[^, which 
is greater than the former estimate and thus satisfies a 
necessary condition on the supply and demand of He. 
Figure 3 shows the GMD of the 13 Gyr isochrones for the 
MS and RGB calculated with two models, one with the 
accretion of He-rich matter ([Fe/H]= -1.3, y=0.24) and 
the other without the accretion ([Fe/H]= -1.6, y=0.24). 
Here, we assume that each star accretes a. Y = 0.4 mat- 
ter at a rate of O.OIA/© Gyr~^ during 1-11 Gyr. For 
reference, stars composed oi aY = 0.4 gas are also plot- 
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Fig. 4. — The relationship between the (B-V) color of ZAHB 
and the envelope mass Mcnv for a fixed He core mass for each 
of the three models (extra mixing model with Vonv = 0.42: solid 
line, Y = 0.24: dashed line, Y = 0.4: dotted line). The He core 
mass for each model is set to be 0.4923, 0.4861, and 0.4632 Mq, 
respectively. In each model, the results with t] = 1/3 in Reimers' 
mass loss formula are included and indicated by arrows. 



Fig. 3. — Theoretical CMDs predicted by our three models, 
i.e., the surface pollution with extra mixing model (solid line), 

Y = 0.24 model (dashed line) and Y = 0.4 model (dotted line) 
together with th e observed CMP for lj Cen, which is indicated by 
the shaded area l|Rev et al.l(200i ). The theoretical CMDs are pre- 
sented by isochrones at the age of 13 Gyr for MS and RGB, and 
evolutionary tracks for HB. The positions of ZAHB are shown by 
various symbols. The extra mixing models with an efficient mass 
loss are denoted by lower two filled circles, attached with evolution- 
ary tracks (solid lines), whereas an open circle is for the model with 
the normal mass loss estimated from t; = 1/3 in Reimers' formula, 
together with an evolutionary track (dot-dashed line). A cross with 
a dotted track is for the Y = 0.4 model, and the other filled circle 
at V~ 15 expresses the Y = 0.24 model. The three ZAHBs in 
the extra mixing model have total masses of 0.61, 0.55, and 0.52 
Mq , respectively, from top to bottom. Here, we adopt the distance 
modulus (m - M)v = 13.5 and the reddening E(B - V) = 0.14. 

ted. It is shown that a. Y = 0.4 matter either on the 
surface (sohd line) or an entire star (dotted hne) shifts 
the color of main-sequence to bluer, as compared to a 
normal Y = 0.24 star (dashed line). 

4.3. Horizontal Branch Morphology 

For later evolution, evolutionary tracks are drawn for 
HB in the CMD from the tip of the RGB. In the sur- 
face pollution model, extra mixing is assumed to operate 
during the RGB and the calculation is performed for a 
O.SMq star, whereas a 0.64 Mq star is considered for the 

Y = 0.4 model because a He-rich star has a short life- 
time. Because of its small mass, the envelope mass of a 

Y = 0.4 HB star is reduced sufficiently to populate an 
extremely blue HB (cross), as dictated by the assump- 
tion o f the standard mass loss formula (77 = 1 /Sl llReimer j 
Il977f) . On the other hand, the zero-age HB (ZAHB) of 
a star, in which the extra mixing operated during the 
RGB, is located less bluer and brighter in the horizontal 
part of HB (open circle), as long as the standard mass 
loss rate is assumed. However, our scenario predicts that 
excess angular momentum in the envelope acquired by 
stellar encounters will enhance the mass loss rate. Thus, 
a star that experienced mixing is expected to have a re- 
duced envelope mass and its predicted location of ZAHB 
becomes close to the observed extreme HB, as shown by 
the lower two filled circles in Figure 3. Such an effi- 
cient mass loss is supported by the existence of dim stars 



aroun d F ~ 19 mag in the observed CMD ijRev et al.l 
l2004f ). Our model predicts that the positions of these 
stars in the CMD are in the pathway of stars that extin- 
guish the H shell-burning during the HB evolution due to 
their decreasing envelopes, and finally move to the white 
dwarf branch. Figure 4 summarizes the dependence of 
the {B — V) color of ZAHB on the envelope mass for 
each model. 

5. CONCLUSIONS 

We first propose that He-rich matter is synthesized in 
red giants with the masses of 0.8 - 2 Mq, which experi- 
ence encounters with other stars in GCs, leading to the 
extra mixing during the RGB phase in their envelopes. 
On the other hand, massive AGB stars are unable to pro- 
duce the He abundance as much as F ~ 0.4 because (1) 
their large envelopes work as a buffer against He enrich- 
ment and (2) the amount of dredged-up matter is limited 
by thinning the He-burning shell, due to the growth of C- 
O core. Therefore, a F ~ 0.4 matter is a product unique 
to GCs. In some GCs, this matter is retained in their 
gravitational potential for a prolonged time and is ac- 
creted by kinematically cool stars, although this may not 
be the case for most GCs. In addition, stellar encounters 
in GCs induce the extra mixing and increase the mass 
loss rate during the RGB phase. As a result, such stars 
are predicted to evolve to significantly blue HB stars. In 
conclusion, MS splitting is inclined to be accompanied 
by the existence of an extremely blue HB, as observed 
in oj Cen and NGC 2808. It should be, however, noted 
that our models might make a variation in the final He 
abundance as a result of extra mixing induced by var- 
ious amounts of angular momentum transfered through 
encounters. If the variation is too large, then it will erase 
the MS splitting and lead to a broad MS. Further inves- 
tigations of the extra mixing model in terms of the MS 
splitting are surely awaited. 

This work is in part supported by a Grant-in-Aid for 
Science Research from the Japansese Society for the Pro- 
motion of Science (19740098). 
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